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EXPONENTIALLY CONVERGENT 
FUNCTIONAL-DISCRETE METHOD FOR 
EIGENVALUE TRANSMISSION 
PROBLEMS WITH DISCONTINUOUS 
FLUX AND POTENTIAL AS A 
FUNCTION IN SPACE L x 




V.L. Makarov • N.O. Rossokhata • D.V. Dragunov 



Abstract — Based on the functional-discrete technique (FD-method), an algorithm for 
eigenvalue transmission problems with discontinuous flux and integrable potential is 
developed. The case of the potential as a function belonging to the functional space L\ 
is studied for both linear and nonlinear eigenvalue problems. The sufficient conditions 
providing superexponential convergence rate of the method were obtained. Numerical 
examples are presented to support the theory. Based on the numerical examples and 
the convergence results, conclusion about analytical properties of eigensolutions for 
nonself-adjoint differential operators is made. 
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1. Introduction 

Physical and chemical processes in multilayer medium can be modeled by partial differential 
equations. Because of different physical and chemical properties of layer materials, besides 
initial and boundary conditions such models involve transmission (matching) conditions 
describing processes in the boundary of two layers. As a rule such conditions provide jump 
of a solution or flux, or both of them. Transmission problem for elliptic equations was 
proposed by M.M. Picone in 1954 [17] for the case of discontinuous solutions and continuous 
flux and then it was studied in [2, 8, 22, 20]. Linear eigenvalue transmission problems were 
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investigated in [3, 4, 16]. The numerical treatment of transmission problems are given in 
[14, 15, 12, 13, 19, 18, 5, 23, 9, 10] et al. Particularly, papers [14, 15, 12, 13, 19, 18] are devoted 
to the development of numerical functional-discrete exponentially convergent methods for 
the eigenvalue transmission problem with discontinuous solutions and continuous flux. 

In this article we extend the FD-approach to the eigenvalue transmission problems with 
continuous solutions and discontinuous flux. 

The paper is organized as follows. Section 2 is devoted to the eigenvalue transmission 
problem with a potential belonging to the functional space Li(0, 1). In section 2 we formulate 
the problem with the additional differential condition providing a uniqueness of the solution, 
then we describe the numerical technique and prove the convergence theorem. In section 3 
the corresponding nonlinear eigenvalue problem has been studied and the convergence result 
similar to that of section 2 has been obtained. The results of numerical experiments are 
discussed in section 4. They are in good agreement with the theoretical findings. Section 5 
includes conclusions about the proposed approach. 

2. Eigenvalue transmission problem with a linear potential as a 
function in L\ 

2.1. Formulation of the problem 

Let us consider the following eigenvalue problem 

u'l(x) + [\-q(x)} Ul (x) = 0, xe(o,^j, q(x)eL 1 (0,l) 

v! 2 \x) + [A - q(x)]u 2 (x) =0, xe Q, , (1) 
with transmission conditions 











= o, 


:••©] 







= 1, 



and boundary conditions 

Ul (0)=«2(l)= 0,^(0) = 1, 
where [/(xo)] = f2( x o) — fi( x o) is a jump of a function / at the point x . 



2.2. Description of the FD-method for the linear case 



According to the FD-approach [11] instead of the original problem (1) we consider the 
following problem with parameter t G [0, 1]: 

d 2 ( 1 

ui(x,t) + [\(t) - tq(x)]u 1 (x,t) = 0, x e 0, - 



dx 2 
dx 2 



u 2 (x, t) + [A(t) - tq(x)]u 2 (x, t) = 0, x G ( -, 1 I , 



9. 



(2) 



«i(0,f) = u 2 (l,t) = 0, t^u'iOM) 



1, 



x=0 
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0, 



1, Vie [0,1], 



where q(x) G Li(0, 1). 

Taking into account that for t = we can easy find the exact solution (A(0), Ui(x, 0), % — 
1,2) to problem (2) and that for t — 1 the solution (X(l),Ui(x, = 1,2) to problem (2) 
coincides with that to problem (1), we can suppose that solution to problem (2) can be found 
in the form of power series with respect to t 



K(t) = 5>£V, u ni (x,t) =J2 u nKx)t j , i = 1,2. 

3=0 3=0 



(3) 



Setting t — 1, we obtain 



A n — A n (l) 



]T u ni (x) = u ni (x, 1) = £ u^x), % = 1, 2 

i=o j=o 



(4) 



provided that these series converge. 

Thus, we can represent the approximate solution to problem (1) as the corresponding 
truncated series 

m m 
An=E A n } ' (X)=£«S ) (X), 1 = 1,2, (5) 

J'=0 j=o 

which are referred to as the approximations of the m-th rank. 

Substituting (3) into (2) and equating the coefficients of equal powers of t, we obtain the 
following recursive sequence of problems: 

for j = — 1 

T(-) + A« ) (x) = 0, *e(o,0, 



u, 



«^(*) + Aff>u$(s) = 0, xef 1 -, 



(6) 



«£ ) (o) = «JS(i) = o,< ; '(o) = i, 



,(o) 



(0)', 



u 



(0) 



0. 



(0)' 



= 1; 



and for j = 0, 1, 2, . . . 

+ Ai%r } (*) = Fif^x, g), x e ( 0, J 



t& +1)W (s) + ^ 2 +1 \x) = F% +1 \x,q), x e Q, 1 



(7) 



«a +i) (o)=«a hi, (i)=o,tff i >)=o, 



,0+i) 



Ci+l)' 



I/, 



Cj+i) 



O'+i)' 



o, 



where 



i^ +1) (a;,g) = - £ A^)^(x) + ^^(z), i = 1,2. 



(8) 



p=0 
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We obviously can find the exact solution (Xn\u^-(x),i = 1, 2), (j = 0, 1, . . .) to each of 
the problems. 

Thus, to apply the proposed algorithm one should find the rough approximation 
(Xn\ Uni(x), i = 1,2) as the exact solution to problem (6) and then consecutively find the 
corrections (Xn\ u^- (x) , i = 1, 2) as the exact solutions to problems (7), (8) for j = 0, 1, 2 . . .. 

2.3. Convergence result 

Formulas (4) and (5) imply that the algorithm's error can be estimated in the following way 



j=m+l 

oo 

||A„-A„||< E 1^1' ( 9 ) 

j=m+l 

provided that series in the right sides are convergent. Therefore, our aim is to find conditions 
providing the convergence of the corresponding series and to estimate their convergence rates. 
To achieve that we will try to construct a geometric sequence with denominator less than 1, 
which dominates the terms of series in formulas (9). 

Firstly, let us find the solution to problem (6). Taking into account the differential 
equation, the boundary and matching, we can write eigenfunction of problem (6) in the 
following form 



sin I \l X^x 



vSHx) = L «£>(*) = cf sin (yfxF\l - x)) , 



where the unknown constant and eigenvalue can be determined using the matching 
conditions 

sm(^X&/2) 

-cf sin ( V / Ai 0) /2) + V == J = 0, (10) 



4 ( Va! — I \/a;,"72) +cosi \/a;;"/2) i. 

From (10) it follows that we have two different sequences of eigensolutions of problem 

(6): 



I. 



Ai 0) =47r 2 (±jj + 2^ , n = 0,l,2... (11) 

V^(l-*)) 



sin Xn^xj sin 

U nli X ) = /= > U n2 \ X ) = 

An V A n 



and 
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II. 



Ai )=4vrV, „£)■ 

'_]\rt+l _ j 

-sin(27rn(l -x)) 



Next, we write the solution to the nonhomogeneous boundary value problem (7) as 



x) 



sin {2Tinx) 
2im 



(12) 



1)™ + 1 

sin (27mx) , n = 1, 2 . . . 



Ci+i) 

nl 



(*) = 



sin 



x(0) 



F,g +1) (Ode, 



(13) 



,Cj+i) 



/ 



i sin (^0* - O) 



(0) 



The unknown constant c^'" 1 " 1 ^ and correction A« of the trial eigenvalue A„ can be found 
from the matching conditions of problem (7). 

Combining representations (13) together with matching conditions (7), we obtain the 
system of linear equations with respect to c^"^ and A„ 



c { i +l) sin 



An 



(0) 



/ v r= — J -F& +1 \t)dt + 

J . U(o) J 

Y An 1/2 



Sill 



,Ci+i) 



cos 



\(0) 
An 



(14) 



/ 



v>s( V^ 0) (l-0 



A 



(o) 



nl 



1/2 



A 



(0) 



n2 



(0« 



Now we are in position to find a recursive formula for the corrections An +1 ^. For this 
purpose we have to consider several cases. 

First, we consider the case when eigensolutions (Xn\u^(x),i = 1,2) are determined by 



formula (11). In this case we have that sin 



v(0) 



7^ and cos 



1 



0. Hence, 
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from system (14) we get the following equation with respect to A„ 



Cj+i) . 



+ 



1/2 

/ 



1 

/ 

1/2 



sin 



COS 



sin 



cos 



sin | \/Ai 0) (|-0 



sin 



cos VWI-O 



cos (y^^ 



d£ = 0, 



which yields us the formula for A« 




1/2 _ 



x 



+ 



1/2 



J2^ +1 - p) uS + q f (x)u^(x) 

p=l 



p=l 



sin ( \J\ { n\l-x) ) dx+ 



sin ( X/X { n\l -x) ) da; 



> . 



Taking into account the equality 



V 2 sin (\f^x) sin (\f^\l - x) ) i sin 2 ( \/Ai 0) (l - x) 

dx+ I - j= y -dx 



(0) 



1/2 



sin I \/An°V2 



2VAf 

we arrive to the following representation 



8 A/ A 



,(o) 



Ai J+1 ) = 



(0) 



1/2 



^Ag +1 -P)«S(x) + 9(x)«S(x) 
p=i 



sin ( \/Al 0) (l -x) ) (15) 



+ 



1/2 



X)Ag +1 ^tig(x) + 9(x)«SW 
p=i 



sin ( \/Ai 0) (l -x) ) cfe 



> . 
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We can find unknown constant c^"^ (see (13)) from either of equations (14). For exam- 
ple, from the first one we have 



sin(VA^/2) 



^sin^V^d-O) 



x(0) 



*fi +1) (0de+ 



- / v p= 

i . A(0) 



1/2 V An 

Formula (16) yields us the estimate 



/ 1/2 



,0+1) 



3Ai 0) 



1 

/ 

1/2 



\ 



i^ +1) (0 



(16) 



(17) 



3Ai 0) 



7 0+i) 

nl 



0,1(0,1/2) 



? 0+i) 

n2 



0,1(1/2,1) 



3A W 



IIpO'+i)!! 

II n 1 1 0,1 ' 



rb 



where || • ||o,i( ,&) denotes the norm in Li(a,6) : ||/(a;)||o,i,(o 1 6) = /„ 1/(^)1^5 and ||/|| 0il = 

fill 0,1(0,1/2) + || /a || 0,1(1/2,1) - 
Hence, formulas (13) - (17) imply the following estimates 



u 



0+i)|| 

n lloo 



fll 0+1) II II 0+1) II 1 ^ 

maxjlKx || oo,(0,l/2), ||<2 ||oo,(l/2,l)j < 



< max 



i r 1/2 



(0) Jo 



\F^ +1 \x)\dx,\^ 



O+i) i 



\F« +1 \x)\dx } ^ 



- -^ll^ +1) (-)|lo,i + ^ll^ +1) (-)llo,i * 



Art 



(0) 



1 ( 2 



2 + >/3 1 

^7! 



i . n u (p) 



n lloo 



+ ll?l|o,i • IK 



0)1 



,p=0 



£|A<? +1 -*> •!!«?) Ik + Iklki • lk j) lloo 



.p=l 



where || • ||oo,M denotes the norm in L^a^b] and WfW^ = max^/iH^o,!^], H/2 ||oo,[i/2,i]}- 
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Thus, to estimate ||t4f' +1 ^||oo and |An~ ri/l |, we have to solve the following system of in- 
equalities 

El A S? +1_P) |- Halloo + Iklkl-ll^lloo 



\u ij+1) \\ < 

\ u n Woo <z 



.P=0 
j 



j'+i-p)| . ||„(p) 



< Hoc + M ^ || 0,1 • ||«. 



0)1 



n lloo 



,p=l 



-i">- ^V A ^ & ^ / ; ;: .. / o.i.2... 



A 



Introducing the new variables 



|u£' +1) |U, and N+l =a>\\^\ 



we can rewrite system (18) in the following form 



(18) 



(19) 



^ ^2vj+i-p v p + IMIo,i u j» 

p=0 



(20) 



Pi+i ^ Vj+i-pVp + IMIo,i u i» 3 = 0,1,2.... 
P =i 

Let us consider the scalar sequences {^}^L > {A*j}>=i defined by the following recurrence 
formulas 



7+1 = 5^j+l-p u P + l|9||o,lUj, 



p=0 

i 

Pj+i = XI Vj+i-pVp + lkl|o,i^j = Uj+i - Jij+iVo, j = 0, 1, 2, 
p=i 

With tJ = t> = ||Mn 0) /^||oo = §■ 

Comparing inequalities (20) with equalities (21), it is easy to see that 

v j+ i ^ V j+ l, fJ, j+ l < JI j+1 . 
Eliminating ~pj + i from system (21) we arrive at the recurrence formula 



(21) 



vj+i = ^2v j+1 - p v p + ||g|| ,iUj(l + uo), j = 0,1,2.... 

P =i 

Multiplying both sides of the last equality by z^ +1 and summing them over j from to 
oo we obtain the following equation 

f(z) -v = [f(z) - v ] 2 + (1 + v ) ||g|| 0) i*/(*) 

or, in more convenient form, 

[f(z) - v } 2 - [1 - (1 + v ) \\q'\\ ,iz] ■ [f(z) - v ] + v (l + v ) \W\\o,iZ = 0, 
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where 



m = E 



Z J Vj. 



(22) 



3=0 



We have the quadratic equation with respect to f(z) — vq with the roots 



(/(*) -Uo),2 = 



where 



L> = (1 +^ ) 2 iklloi x 



(23) 



x 



1 + 2^ + 2^0(1+^0) \ / 1 + 2u - 2^0(1+^0) 



(1+Uo) lkllo,l 



v! + Uo) lkll ,i 



It is easy to see that the solution which represents the generating function (22) is the following 



one 



/ (z) = v + 



(l-(l+77 ) IMIo,l^) "V^ 



(24) 



It is obvious that the right-hand side of equality (24) can be expanded as a power series in 
z where z G [0, R] and 



R = 



1 + 2^0-2^0(1 + ^0) 



' X + H ^ II 0,1 (1 + v ) \\q\\ 01 (l + 2v + 2y/v (l + v j) 



(1 + 1) IMIo.i + f + 



(25) 



Here i? denotes the radius of convergence for series (22). On the other hand, it follows from 
convergence of series (22) that 

R j v j 



for some positive constants c and e. 

Returning to (19), we arrive at the estimates 



u 



(i+i) 1 



3+1 



3c 



8(J + 1)^V^ 



(26) 



and 



i\(i+i)i = Mj+i < /fj+i < gj+i < c 

ln 1 a* ^ ^ (1 + uoV ^ (j + l) 1+£ R 



V3VK 
2+^3 



(0) 



R 



(27) 



Thus, series (4) converge, and the last two inequalities imply the error estimates 

m+l 



I _ m i < II 0')|| <r 

1 1 00 ^; / || M n ||oo 
j=m+l 



3C 



2 + \/3 



8(7^ + 1)^7^ \V3\f^R 
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c 



|A„- A„ | < Y. 1^1 < taTI 



j=m+l 



provided that 



2 + V ~ 3 <i, 



2 + \/3 



where 



C 



2+V3 



(28) 



(29) 



Substituting expression for i? (25) into formula (29) we arrive at the following convergence 
condition for the algorithm (5), (6), (7) 



(2 + V3)£||?llo,i (l + f + 2y/f) 
^= w -^-<l. 



,(0) 



(30) 



Now let us consider the case when eigensolutions (Xn\u^(x),i = 1,2) are determined 
by formulas (12). We see that in this case sin A^/2^ = 0. Hence, to find An +1 ^ we 

should use the first equation of (14) and parameter c$ +1 ^ can be determined from the second 
equation of (14), that is, 



1/2 1 

A 0'+i) _ ) f sin2 ( 2nnx )^ i ii 1 mi i -\\ f sm 2 (27mx) 



2irn 



-dx + ((-l) n + 1) 



2-nn 



dx > x 



(31) 



1/2 



1/2 



X 



J2^ +1 - p) u^(x) + q (x)uS(x) 



p=l 



sin {2imx) dx+ 



+ 



1/2 



-^AS" +1 -rt«g(x) + 9(x)«g(x) 



p=i 



sin (27rna;) cte > , 



U+n _ /" 1/2 cos (2™) Z" 1 cos(27rnx) 0+1) 



j+i) _ _ f 
1 ~ Jo 



2nn 



-F^ +l, (x)dx 



1/2 



2-KTL 



Fg^\x)dx. (32) 



It is easy to verify that 



1/2 

r- 



sin 2 (2nnx) 



2ixn 



dx + ((-l) n + l) 



sin 2 (2nnx) 
2-nn 



dx = < 



1/2 



8"7m 
1 

8™ 



, n = 2k, k G N, 
, n = 2k + 1, k E N. 
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Co-, sin (2nnx) 

If n is an even number, then u n2 :(x) = . Substituting the last expression into 



formula (32), we get 



irn 



8im 



1/2 



sin (2nnx) dx+ (33) 



+ 



1/2 



p=l 



sin (2irnx) dx 



Formulas (13), (32) and (33) imply that 



l4 J+1) Hoo = max | ||oo,(0,1/2), l|oo,(l/2,l) } < 



and 



< max 



-r= r\F^\ x) \d X ,\cr ) \ + ^= f |fg +1 >. 



. 1 1 

^ I + 



■Inn 2tt/; /' 1 "'^ (x) Ho,i ^ 



7T71 



El^^l-ll^lloo +||?||0,l||«S 



00 1 



p=0 



|A^' +1) | < 



87m 



EM 



i . \\ u (p) 



n lloo 



+ IMkilK J 



that is, we get inequalities (18) with 



(34) 



(35) 



a = nn, b = 



8irn 



Repeating the same computations as for the previous case, we obtain the similar error 
estimates 



provided that 



where 



K^n ^rt ||oo ^ 



3C 



m 

An.— A„ ^ 



8nn(m + l) 1+e \7rnR 
C / 1 \ 



m+l 



r n = 



n 

1 T 
irnR 



(m + l) 1+£ R \7inRJ 
kilo,! (l + f + 2^) 

^ v — ^ < 1, 



C 



(36) 

(37) 
(38) 



14 



V.L. Makarov et al. 



If n is an odd number, then inequality (35) should be replaced by the following one 



|A^ +1) | < 8nn 



ElAg+^l-II^Hoo + IMkllltt 

.p=l 



0) II 

n lloo 



(39) 



however, the estimate for ||ttn +1 ^||oo remains the same, see (34). Thus we arrive at the 
equation (18) with 

a = Tin, b = 

and we get the following error estimates 



C / 1 \ m+1 

8iin(m + l) 1+e \7mR J 



U n || OO ^ o | I r, ": 



m C ( 1 



with C defined in (38), provided that inequality (37) holds true. 
Therefore, we have obtained the following convergence result. 

Theorem 2.1. 



(a) Assume that the index n of a trial eigenpair satisfies condition (30) with A„ 
47T 2 (±| + 2n) 2 ; {n = 0, 1, . . .). Then the numerical algorithm (5), (11), (13), (8), (15), 
(16) converges to the corresponding eigensolution of problem (1) super- exponentially 
with error estimate (28). 

(b) Assume that the index n of a trial eigenpair satisfies condition (37) with = 4ir 2 n 2 , 
n = 1,2,. . . . Then the numerical algorithm (5), (12), (13), (8), (31), (32) converges to 
the corresponding eigensolution of problem (1) super- exponentially with error estimate 
(36) for even number n and with error estimate (40) for odd number n. 

Using the term "super-exponential convergence", we intend to emphasize the fact that 
owing to the presence of factor (m + l) 1+e in the denominators of estimates (28), (36) and 
(40) the method remains convergent even if r n — 1. 

3. Eigenvalue transmission problem with a nonlinear potential in- 
cluding a function in space L\ 

3.1. Formulation of the problem 

Let us consider the following eigenvalue problem 

u'l(x) + [A - q(x)] Ul (x) - N( Ul (x)) = 0, xe(o,^j, q(x) & £i(0, 1) 

u'^x) + [A - q{x)]u 2 {x) - N(u 2 (x)) =0, x G Q, , (41) 
ui(0) = « 2 (l) = 0,«i(0) = l, 
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u 



0, 



u 



(42) 





d n 


, N{u) 




du n v ' 


du n v 7 



where A" : R — > R is an analytical function with respect to u, N(0) = 0, that is, 

oo 

N(u) = a iu\ Vm el, di G R, Vi G N. 
i=i 

By A" (w) we denote a fixed function from C°° (R) such that 

, n = 0,l,.... 

«=|«| 

Such a function exists. Moreover, as it follows from (42), the function N (u) 
satisfies inequalities (43) for any n e N. 

3.2. Description of the FD-method for the nonlinear case 

Using the general idea of FD-method, we approximate a numerical solution to problem 
(41) by the truncated series (5) where an eigenpair (Xn\u^(x),i = 1,2) is the solution to 
problem (6) (see formulas (11) and (12)) and pairs u ( i +1 \x), i = 1, 2) (J = 0, 1,2,.. .) 

are the solutions to the nonhomogeneous transmission problems (7) with 

F^ 1] (x, q) = -J2 ^ +1 - P) u { :}{x) + q{x)u%{x) + A 3 {N-u^ , ...,«£>), (44) 

p=0 



(43) 



oo 

i=i 



where 



(45) 



E ^ Ql) (^ } (^)) 



,0-1) 



ai+...+«j=j 
ai>...^aj>0 



(ai - a 2 )\ 



(a,-)! 



, J>0, 



= ^(«2 ) ) 

are the Adomian polynomials (see [21], [1]). 
3.3. Convergence result 

As it was in the linear case described in section 2, zero approximation (Xn\u^(x),i = 1,2) 
can be found using formulas (11) and (12). Applying the same approach as in section 2, we 
determine u~^ l \x) (i = 1, 2) according to formulas (13), (44), (45). The unknown constants 
4 J+1) and can be found using the matching conditions (14), (44), (45). 

Hence, for the case of eigenpairs (Xn\u^(x),i = 1,2) determined by formulas (11), we 
obtain the following expression for Xn +1 ^ 



(46) 



16 



V.L. Makarov et al. 



.1/2 



X 



p=i 



nl y 



sin a/ An°^(l — :r)Gfe+ 



+ 



1/2 



P =i 



sin v/ A^ - 1 (1 — V , 



and parameter c^ 4 ^ is determined by formula (16) with (44). 

Similarly to the linear case discussed in section 2, formulas (13), (16), (44) together with 
(42) and (45) yield us the following system of inequalities 



u 



Cj'+i)ii < 

n Woo ^ 



EM 

,p=0 



j+i-p)| . || u (p) 



n Moo 



+ ll?lkl-||tig ) ||oo+ 



+ E ^ ai) (ik o) iioo) 



(!) I 



«1— «2 



11^-1—^ |U,0')|| Q i 

«n oo oo 



a\+...+aj=i 
a 1 ^...^a j ^0 
a;6Z 



(ati - a 2 )\ 



(aj_i-Q!j)! («j)! 



(47) 



+ E ^ ai) (ii 

ai+...+aj=j 



Ei^ +1_p) i-ii^ ) iioo+ikiki-ik j) iioo+ 

,p=l 



(0)11 A 



L,(l) llai— c*2 
l u ™ Moo 



!) M^j-i— |U,0')|| Q i 
"n oo u-n oo 



(ai-a 2 )! (aj_i-aj) ! 



with a 



V3_J\(P) h - 3 1 



Ar, & 



2+v^ » u ^ A 

Introducing new variables 



(0) 



^■+i = -i-||«{? +1) lloo and ^ +1 = o»-|Ag +1 )|, J>-1, 



(48) 



we arrive at the following system of inequalities: 



< ^2h+i-p v p+ hhivj + A j {N;v ,v 1 ,...,v j ), 

p=0 



(49) 



^'+! ^ E^' +1 -p^ + II9||o,iWj + ^4j-(iV; v , vi, --, 3 = 1,2,.... 



p=i 



To obtain similar estimates for t> 1 and ji\ let us consider inequalities (47) with j = in more 
detail: 

' "'^l + lklkOll^L+^dhfllJ 



.mi 



|W n ' ||oo ^ 



ia^Ioo^MnIo^KIL+^ckiL) 



(50) 
(51) 
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Using notations (48) from (50) we obtain 

Vl < (/^l + lkllo,l) U 0+^iV(||4 0) ||oo) = 

= (/ii + lkllo,i)^o + ^(iV(ll4 0) lloo)-iV(0))= (52) 

= + Ikllo.i) v ° + \n' [ewufw^) iu < 

< + IMIo,i) Vo + N'(v )v . 
Similarly to (52) we can get the following inequality for 

A*i < || ^ || o,i v o + N' (v ) v . (53) 

Now let us consider the sequences {vj}^ and {/Z -}?1 defined by the following recurrence 
equalities 

3 

Vj+i = ^Jl j+1 - p Vp+ \\q\\o,iVj + Aj(N;v ,vi, ...,Vj), (54) 
j 

JI j+l = ^Ji j+l _ p Vp + ||g|| 0) iUj + Aj( N '^o,vi, ...,Vj) = v j+1 - JI j+1 vo (55) 
P =i 

for j'gN and 

«i = (a*i + Ikllo.i) + N' (v ) v , (56) 

A*i = IMIo,i + ^' (vo) v = V!- JijVo, (57) 

where vo = vo = §. It is easy to see that Vj+i ^ ?7j+i and ^ V j G N U {0} . 

Eliminating variables from equalities (54)-(57) we arrive at the recurrence formulas 

3 

vj+i = ^ v j+1 -pV p + (l + vo) (||g||o,iUj + Aj(N;v ,vi, ...,Vj)) (58) 
P =i 

for j G N and 

vi = (1 + uo) (lkllo,i + A 7 ' (v ) vo) ■ (59) 

Applying the technique similar to that used in section 2 to equalities (58), (59) we can 
obtain the following nonlinear equation with respect to generating function f(z) 

f(z) -v = [f(z) - v ] 2 + (1 + v )z[\\q\\ 0>1 f(z) + N(f(z)) + N' (v ) v -N (v )}, (60) 

oo 

/(*) - (61) 

j=0 

Let us prove that power series (61) possesses the nonzero radius of convergence. For this 
purpose we consider the inverse mapping z = f~ l . From equation (60) it follows that 

z (f) = f-^-lf-Jo? = (62) 

(l + v )[\\q\\o,if + N(f) + N (vo)v -N(v )} 
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Since function z(f) (62) is holomorphic in some open interval containing the point fo = Vq 
and z (vq) = 0, we can easily calculate the value of derivative z'(vq) 



it \ ,. z(f)-z(y ) 
z (v ) = hm - — 

f->v J - Vq 

= lim 



(63) 



l-(f-vo) 



(1 + uo)[||g||o,i/ + N(f) + N (v ) v -N (v Q )] 
1 

> 0. 



uo(l + uo)[||g||o,i + ^ (*><>)] 

Inequality (63) implies that there exists an inverse function f — z -1 which is holomorphic 
in some interval (—R,R) (see [7, p. 87]). Now let us prove that series (61) converges at 
the endpoint z = R. Suppose that this is not the case and series (61) diverges at the point 
z = R, that is, 

lim f(z) = +oo. 

However, taking into account that equality (60) holds for every z in (-R, R), we immediately 
get the contradiction 

1= lim (f(z)-v + 

z— >R— 



(1 + Vo)z[\\q\\ 0>1 f(z) + N(f(z)) + N (v ) v - N (go) ] 

f(z) - Vq 



+oo. 



(64) 



This contradiction implies that f(R) < +oo. Thus, we have the inequality 

R j v j < 



;l+e 



J 

with some positive constants c and e, where the constant R is determined by the value of 
||?||o,i and function N (u) . 

Returning to notation (48), we formally obtain estimates (26) and (27). Hence, for n 
sufficiently large the numerical solution converges to the exact one with error estimates (28), 
where R denotes the radius of convergence to the power series (60), (61). 

Now we apply the same technique as above to the case when Xn^ = 4n 2 n 2 (n = 1, 2, . . .). 

If n is an even number then 



8nn 



1/2 



P =i 



sin (2nnx) dx— (65) 



1/2 



P =i 



sin (2imx)dx > . 



and functions u^}(x) (i = 1, 2) can be found using formulas (13), (44) where parameter c% 
is determined by formulas (32), (44). 
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If n is an odd number, then 

A 0'+i) = 

V2 r , 



8nn 



P =i 



sin {2imx) dx— (66) 



1/2 



^2 X { J +1 ~ p) u^ + q(x)u { H(x) + Aj(N; u { ° 2 \ ...,u 
P =i 



nil 



sin (2imx)dx 



and functions u^-(x) (i = 1, 2) can be found using formulas (13), (44) where parameter c. 
is determined by formulas (32), (44). 

As a result we formally obtain estimates (36) for an even number n, and estimates (40) 
for an odd number n, provided that 



Cj+i) 



< 1. 



irnR 

Therefore, we arrive at the following convergence result for the nonlinear problem. 
Theorem 3.1. 



(67) 



(a) Suppose that the index n of a trial eigenpair satisfies condition (29) where = 
47T 2 (±| + 2n) 2 , (n = 0, 1, . . .) and R is a radius of convergence to series (61) satisfying 
equation (60). Then FD-algorithm (5), (11), (44), (46) converges to the corresponding 
eigensolution of problem (41) super- exponentially with error estimate (28). 

(b) Suppose that the index n of a trial eigenpair satisfies condition (67) where R ia a radius 
of convergence to series (61) satisfying equation (60) and A^ = 4ir 2 n 2 (n = 1,2, . . .). 
Then the FD-algorithm (5), (12), (44) an d (65) (if n is even) or (66) (if n is odd) 
converges to the corresponding eigensolution of problem (41) super- exponentially with 
error estimate (36) for even n and with error estimate (40) for odd n. 

4. Numerical examples 
4.1. Example 1 

As an example we consider the following nonlinear eigenvalue problem 

u"(x) -(x + 3x 2 )mi(x) + Xu!(x) - [u^x]} 2 = 0, x e ^0, 

u 2 \x) — (x + 3x 2 )u 2 (x) + Xu 2 (x) - [u 2 (x)] 2 = 0, x E l^j 



(68) 



u 



0, 



u 



ui(0) =u 2 (l) = 0, u'(0) = l. 
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Let us compute the approximations to the six least eigenvalues of the problem. From 
formulas (11) and (12) we obtain the zero approximations for the FD-method's algorithm: 
from (11) with n = : 



A 



(0) 



16 
~9~ 



2 (0) / \ 

7T , u"{ (x) 



Mfa) ro) ( , 

4 ' "0,2 \ X J 



sin (|tt (1 - x)) 



from (12) with n = 1 : 
from (11) with n — 1 : 



sin(27ra;) 



2tt 



$ (s) = 0; 



x(0) 
A 2 



64 
~9~ 



-7T 



sin(|7ra;) 



«2?2 fa) 



7T 



sin(|7r(l-x)) 



from (12) with n = 2 : 

A<°> = l&r 2 , 4°i» = 
from (11) with n — 1 : 

«fi 0*0 

for (12) with n = 3 : 



256 2 



sin(47rx) 
An 



sin(y7ra;) 



u 



sin (4tt (x)) 
4tt ' 



16 



U 



4,2 ( X ) — 



sin (f 7T (1 - a;)) 



16 



7T 



A?) = 36^, = !^, «S(x) = 0. 



(69) 



(70) 



(71) 



(72) 



(73) 



(74) 



We have calculated the FD-approximations for the eigenvalues A n and corresponding 
eigenfunctions u n (x), n = 0, 5 up to the rank m = 4 inclusive. Since the exact solution of 
problem (68) is unknown, in order to estimate the error, we analyzed the Loo-norm of the 
residual 

Vni{x) = -j-5 u m{x) - {x + 3x )u ni (x) + \ n u ni {x) - |u ni (x)J . (75) 
The results obtained are presented in tables 1-7 and depicted in figures 1-5. 



Table 1. Example 4.1. FD-approximation of A and uq(x). 



m 


m 

A 


ll M !u(>)IU[0,l/2] 


IkffWIloo,^,!] 


||^ m) (a:)||oo,[o,i] 





17.5459633797144 


0.24 


0.24 


0.55 


1 


19.6940699073641 


0.19e-l 


0.25e-l 


0.19e-l 


2 


19.6740808021547 


0.60e-3 


0.11e-2 


0.96e-3 


3 


19.6754846046439 


0.28e-4 


0.51e-4 


0.36e-4 


4 


19.6754786167117 


0.88e-6 


0.20e-5 


0.22e-5 



1 



A = 19.6754786167117. 
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x 

a) b) 
Figure 1. Example 4.1. FD-approximation of uq(x) (a) and 1(2(2;) (b). 



Table 2. Example 4.1. FD-approximation of Ai and u\(x). 



m 


m 

Ai 


OOIU,[0,l/2] 


Ikff^lU.tiAi] 


lki m) (^)lloo,[0,l] 





39.4784176043574 


0.16 


0. 


0.11 


1 


40.0755170720146 


0.80e-3 


0.13e-2 


0.26e-2 


2 


40.0597952320099 


0.94e-4 


0.93e-4 


0.93e-4 


3 


40.0595106843757 


0.20e-5 


0.29e-5 


0.63e-5 


4 


40.0594734299829 


0.27e-6 


0.27e-6 


0.42e-6 



4 



Ai = 40.0594734299829. 




a) b) 
Figure 2. Example 4.1 FD-approximation of U\(x) in [0, 1] (a) and in (0.5, 1] (b) 
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Table 3. Example 4.1. FD-approximation of A 2 and U2(x). 



m 


TO 

A 2 


IK™ ( a; )lloo,[0,l/2] 


lk£?( x )IU[i/2,i] 


\\ v 2 Op)||oo,[0,1] 





70.1838535188575 


0.12 


0.12 


0.36 


1 


71.9666624409054 


0.37e-2 


0.46e-2 


0.48e-2 


2 


71.9768253657124 


0.53e-5 


0.86e-4 


0.67e-4 


3 


71.9766663564735 


0.12e-5 


0.20e-5 


0.23e-5 


4 


71.9766690902938 


0.42e-7 


0.46e-7 


0.41e-7 



1 



A 2 = 71.9766690902938. 




a) b) 
Figure 3. Example 4.1. FD-approximation of 1*3(0;) (a)) and Ui{x) (b)). 



Table 4. Example 4.1. FD-approximation of A3 and u${x). 



m 


TO 
A 3 


ll«3a O c )lloo I [0,l/2] 


Ikff^lloo.iiAi] 


||^3 ro) (^)||oo,[o,i] 





157.913670417429 


0.80e-l 


0.16 


0.49 


1 


159.737504889796 


0.13e-2 


0.20e-2 


0.17e-2 


2 


159.739358889030 


0.18e-4 


0.18e-4 


0.22e-4 


3 


159.739350000888 


0.22e-6 


0.21e-6 


0.21e-6 


4 


159.739350058922 


0.48e-9 


0.85e-9 


0.15e-8 



4 



A 3 = 159.739350058922. 



Table 5. Example 4.1. FD-approximation of A4 and 114(2:). 



m 


TO 

A 4 


||«4a ) (aj)lloo,[o I i/2] 


IkS ( z )IU.[V2,i] 


||^ m) (^)||oo,[0,l] 





280.735414075430 


0.60e-l 


0.60e-l 


0.21 


1 


282.620725160874 


0.57e-3 


0.13e-2 


0.18e-2 


2 


282.622515077477 


0.12e-4 


0.10e-4 


0.11e-4 


3 


282.622528329338 


0.36e-7 


0.79e-7 


0.56e-7 


4 


282.622528620046 


0.28e-9 


0.24e-9 


0.59e-9 



i 



A 4 = 282.622528620046. 
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Table 6. Example 4.1 FD-approximation of A5 and 115(2:). 



m 


m 

A 5 


||«5?( x )lloo,[o,i/2] 


||4?(>)IU,[i/2,i] 


lM m) (>)l|oo,[o,i] 





355.305758439216 


0.53e-l 


0. 


0.53e-l 


1 


355.816547268956 


0.11e-3 


0.47e-4 


0.14e-3 


2 


355.814884402534 


0.83e-6 


0.10e-5 


0.15e-5 


3 


355.814878976396 


0.24e-8 


0.96e-8 


0.47e-7 


4 


355.814878544097 


0.20e-9 


0.30e-9 


0.59e-9 



4 



A 5 = 355.814878544097. 




a) b) 
Figure 4. Example 4.1. FD-approximation of 1*5(2;) in [0,1] (a) and in (0.5,1] (b). 
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Analyzing the results of calculations we have noticed that there exist eigenfunctions with 

some indexes n such that their numerical approximation u n (x) have exactly n — 1 zeros in 
the interval (0,1) (see figure 2 for n = 1 and figure 3 for n = 5). Zero approximations 
of such eigenpairs are determined by formula (12) with odd n. Thus, taking into account 
the numerical calculations and convergence result, we can conclude that there does not 
exist an eigenfunction of problem (68) with exactly one or five zeros in the interval (0, 1). 
Further numerical experiments suggests that there does not exist an eigenfunction which has 
exactly 1 + 4k (k — 0, 1, 2, . . .) zeros in the interval (0, 1). Therefore, there exist two linearly 
independent eigenfunctions with exactly 4k zeros. 

Hence, we can conclude that the theoretical approach developed in [6] and [24] for non- 
linear eigenvalue problems is not applicable for the problems of type (1). The reason is that 
the corresponding differential operator is nonself-adjoint. So we see the crucial importance of 
the requirement that the differential operator is self-adjoint for the applicability of Theorem 
4.2 and Theorem 5.6 from [6]. 

To illustrate the exponential convergence rate of the method it is convenient to consider 
the function 

I = l n (m) = ln( Halloo). (76) 

The values presented in table 7 suggests that the graphs of function (76) for n G 0,5 are very 
similar to the straight lines, see Figure 5. And it is this fact that confirms the exponential 
nature of the FD-method's convergence. Also, it is easy to see that the slope of the lines 
on Figure 5 increases as the index n of the trial eigenvalue increases. This means that the 
convergence rate of the method increases along with the index of eigenvalue. 

Table 7. Example 4.1. The values of function l n (m) (76) 





n = 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


171 = 


-0.60 


-2.21 


-1.02 


-0.71 


-1.56 


-2.94 


m = 1 


-3.96 


-5.95 


-5.34 


-6.38 


-6.32 


-8.87 


m = 2 


-6.95 


-9.28 


-9.61 


-10.72 


-11.42 


-13.41 


m = 3 


-10.23 


-11.97 


-12.98 


-15.38 


-16.70 


-16.87 


m = 4 


-13.03 


-14.68 


-17.01 


-20.32 


-21.25 


-21.25 



4.2. Example 2 

As the second example let us consider the following nonlinear eigenvalue transmission prob- 
lem with singularity 

u'[{x) + (x - {\ - x)- 1 ^ Ul (x) - M*)] 2 = 0, x G (o, ^ , 

v»(x) +(\-(x- l)- 1 ^ u 2 (x) - [u 2 (x)] 2 = 0, x G Q, l) , (77) 
ui(0) = « 2 (1) = 0, u' 1 (0) = 0. 
It is easy to ensure that for the case when q(x) = — -== the indefinite integrals in 



2 ^1 



formulas (13) cannot be expressed through the elementary functions. Therefore, to approx- 
imate functions u^{x) we have used the numerical integration according to the Simpson's 
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Figure 5. Example 4.1. Loo-norm of residual (75) versus the rank m of the approximation on a semi- 
logarithmic scale. 

rule with the precision 10 -15 (for n = 0, 1, 2) and 10~ 19 (for n = 3) 1 . Since the exact solution 

to the problem is unknown and u n i (|) = oo, for the error control we have used the norm of 
residual 



/Ml 



,( m )i 



,0)i 



W n ||oo(o,i) — m ax <j ||z/ n)1 || 00(0,1), lloo(|,i) 



where 



dx 2 



Unlix) + \ r 



— X 



(7* 



d 

dx 



Unl(x) + 



X 


[( 


I 







An 



ra , . ,m , . . . 
Mnl(^) - {u nl (x)Y 



X 



dx, x e (0, 



d 

U. m 



dx A 



Un2(x) + X n 



m , . / m / 
Wn2W - {Un2{X)) 



X 



dx 



(m) 



x e 



1 \ dm dm /1\ 

2 + ^" 2(x) -^ 2( 2 ) + 



A,, 



U n 2{X) ~ {U n2 {x)y 



X 



dx, 



The results of calculations for the least four eigenvalues and corresponding eigenfunctions 
are presented in Tables 8-11 and Figures 6-9. 

Similarly to Example 4.1, analyzing the results of calculations we have observed that 
there also exist eigenfunctions with some indexes n such that their numerical approximation 



The precision was estimated using a posteriori error analysis. 
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Table 8. Example 4.2. FD-approximation of A and Uq(x). 



171 


m 


II ( m )/ Ml 

\ \ H 1 1 Of* 1 r „ _. i „n 

II M 0,1 WI|oo,[0,l/2] 


II ( m ) 1 Ml 

|P0,2 WI|oo,[l/2,l] 


ii ( m ) ii 

11^0 loo, [0,1] 


n 

w 


1 7 54^96337971 44D1 


24 


24 


^8 


l 


21.814604661812502 


0.24e-l 


0.24e-l 


0.91e-2 


2 


21.733545015000821 


0.15e-2 


0.15e-2 


0.33e-3 


3 


21.734895246851330 


0.67e-4 


0.67e-4 


0.16e-4 


4 


21.734885594786305 


0.24e-5 


0.24e-5 


0.49e-6 


5 


21.734887489687971 


0.90e-7 


0.90e-7 


0.22e-7 


6 


21.734887360755933 


0.42e-8 


0.42e-8 


0.11e-8 


7 


21.734887362837029 


0.18e-9 


0.18e-9 


0.51e-10 


8 


21.734887362829545 


0.73e-ll 


0.73e-ll 


0.18e-ll 



8 



A = 21.734887362829545. 



Table 9. Example 4.2. FD-approximation of Ai and u\(x). 



m 


m 

Ai 


ll w u( a Olloo,[0 ) l/2] 


lki?(^)||oo,[i/2,i] 


ii Mn 

\\ V 1 Hoc, [0,1] 





39.478417604357434 


0.16 





0.12 


1 


41.751445051880136 


0.15e-2 


0.31e-2 


0.28e-3 


2 


41.751095410924888 


0.17e-4 


0.26e-4 


0.13e-4 


3 


41.751103408518688 


0.69e-7 


0.20e-5 


0.19e-6 


4 


41.751101792373055 


0.10e-7 


0.30e-7 


0.91e-8 


5 


41.751101774628877 


0.19e-9 


0.12e-8 


0.23e-9 


6 


41.751101775581723 


0.63e-ll 


0.35e-10 


0.59e-ll 


7 


41.751101775606647 


0.21e-12 


0.62e-12 


0.21e-12 


8 


41.751101775606187 


0.33e-14 


0.33e-13 


0.33e-14 



\i = 41.751101775606187. 



Table 10. Example 4.2. FD-approximation of A2 and u 2 (x). 



m 


m 

A 2 


1 1 Uj? ( X ) II 00, [0,1/2] 


IkffWIIooJiAl] 


11 Mn 

\\ U 2 Hoc, [0,1] 





70.183853518857661 


0.12 


0.12 


0.55e-l 


1 


73.434236059947037 


0.36e-2 


0.36e-2 


0.57e-3 


2 


73.462388394109106 


0.39e-4 


0.39e-4 


0.78e-4 


3 


73.462228765029615 


0.33e-5 


0.33e-5 


0.18e-5 


4 


73.462193970765803 


0.92e-7 


0.92e-7 


0.11e-6 


5 


73.462193827166322 


0.52e-8 


0.52e-8 


0.48e-8 


6 


73.462193886381801 


0.21e-9 


0.21e-9 


0.16e-9 


7 


73.462193887216886 


0.89e-ll 


0.89e-ll 


0.13e-10 


8 


73.462193887097591 


0.54e-12 


0.54e-12 


0.24e-12 



A § 2 = 73.462193887097591. 
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Table 11. Example 4.2. FD-approximation of A3 and 113 (x). 



m 


m 
A 3 


ll w 34 0)||oo,[0,l/2] 


ll u 3?( a: )lloo,[l/2,l] 


Ir3 Hoc, [0,1] 





157.9136704174297379014 


0.89e-l 


0.16 


0.66e-l 


1 


160.2464778440802534274 


0.50e-3 


0.14e-2 


0.88e-4 


2 


160.2466449084238753618 


0.18e-5 


0.26e-5 


0.13e-5 


3 


160.2466411888247570166 


0.33e-8 


0.62e-7 


0.17e-7 


4 


160.2466412105927434146 


0.47e-10 


0.63e-9 


0.11e-9 


5 


160.2466412109016291827 


0.50e-12 


0.33e-ll 


0.12e-ll 


6 


160.2466412109020454979 


0.33e-14 


0.37e-13 


0.12e-13 


7 


160.2466412109020583699 


0.31e-16 


0.42e-15 


0.11e-15 


8 


160.2466412109020585522 


0.36e-18 


0.34e-17 


0.33e-17 



8 



A 3 = 160.2466412109020585522. 





a) b) 
Figure 7. Example 4.2. FD-approximation of u\{x) in [0,1] (a) and U\{x) in (0.5,1] (b). 
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Figure 8. Example 4.2. FD-approximation of 113(0;). 




Figure 9. Example 4.2. Loo-norm of Vn (x) versus the rank m of the approximation on a semi-logarithmic 
scale. ^: n = 0, •: n = 1, ■ : n = 2, o: n = 3. 
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u n (x) have exactly n+l zeros in the interval (0, 1) (see figure 7 for n = 1). As in Example 4.1, 
zero approximations of such eigenpairs are determined by formula (12) with odd n. Thus, 
taking into account the numerical calculations and convergence result, we can conclude that 
there does not exist an eigenfunction of problem (69) with exactly Ak + 1 (k — 0, 1, 2, . . .) 
zeros in the interval (0,1). Moreover, there exist two linearly independent eigenfunctions 
with exactly Ak + 2 zeros. 

Thus, Example 4.2 also demonstrates un-applicability of analytical theory developed in 
[6, 24] for the case of nonself- adjoint differential operator. 

The numerical experiment also confirms the exponential convergence rate of the proposed 
method improving along with the increase of the index of a trial eigenvalue (see figure 9). 



5. Conclusions 



In the paper we present a recursive algorithm for solving both linear and nonlinear eigen- 
value transmission problems with potential belonging to space Li(0, 1), discontinuous flux 
and continuous solution. Using the method of generating functions, we find the sufficient 
conditions providing the exponential convergence rate of the method. Also, we establish 
the fact that the convergence rate of the method improves along with the increase of the 
index of a trial eigenpair. Numerical results confirm the theoretical conclusions. Further- 
more, numerical examples and convergence results demonstrate the crucial importance of 
the assumption that differential operator is self-adjoint for the applicability of the theory for 
nonlinear eigenvalue problems developed in [6] and [24]. As it follows from the numerical 
examples, the fact that differential operator is nonself- adjoint results in the lack of strong 
dependence between the index of an eigenpair and the numbers of zeros of the corresponding 
eigenfunction. 

References 

[1] K. Abbaoui, Y. Cherruault, and V. Seng, Practical formulae for the calculus of multivariable adomian 
polynomials, Math. Comput. Modelling, 22 (1995), no. 1, pp. 89-93. 

[2] S. Campanato, Sui problemi al contorno per sistemi di cquazioni differentiali lineari del tipo 
dell'elasticita, Annali della Scuola Normale Superiore di Pisa: Scienze Fisiche e matematiche, 13 (1959), 
no. 2, pp. 223-258. 

[3] M. N. Ozisik, Boundary value problems of heat conduction, International Textbook Company, Scranton, 
Pensylvania, 1968, 506P. 

[4] M. N. Ozisik, Heat conduction, John Wiley and Sons, New York, 1980, 688P. 

[5] D. Givoli, Finite element modeling of thin layers, Computational Modeling of Engineering Sciences, 5 
(2004), no. 6, pp. 497-514. 

[6] H. Heinz, Nodal properties and variational characterizations of solutions to nonlinear Sturm- Liouville 
problems, Journal of Differential Equations, 62 (1986), pp. 299-333. 

[7] E. Hille, Analytic function theory, Ginn and Co., Boston, 1959, 308pp. 

[8] L. J.L., Contributions a un probliemc de M.M. Picone, Annali di Matematika Pura ed Applicata, 41 
(1955), no. 4, pp. 201-219. 

[9] B. Jovanovich and L. Vulkov, Finite difference approximation of strong solutions of a parabolic interface 
problem on disconnected domains, Publications of Institute of Mathematics, 84 (2008). 

[10] B. Jovanovich and L. Vulkov, Numerical solution of a hyperbolic transmission problem, Computational 
Methods in Applied Mathematics, 8 (2008). 



30 



V.L. Makarov et al. 



[11] V. Makarov, About functional-discrete method of arbitrary accuracy order for solving sturm-liouville 
problem with piecewise smooth coefficients, Dokl. Akad. Nauk. SSSR, 320 (1991), no. 1, pp. 34-39. 

[12] V. Makarov and N. Rossokhata, FD-method for nonlinear eigenvalue problems with discontinuous eigen- 
functions, Nonlinear Oscillations, 10 (2007), no. 1, pp. 126-143. 

[13] V. Makarov and N. Rossokhata, A rewiev of functional-discrete technique for eigenvalue problems, 
Journal of Numerical and Applied Mathematics, 97 (2009), pp. 97-102. 

[14] V. Makarov, N. Rossokhata, and B. Bandyrskii, Functional- discrete method with a high order of accuracy 
for the eigenvalue transmission problem, Computational Methods in Applied Mathematics, 4 (2004), 
no. 3, pp. 369-381. 

[15] V. Makarov, N. Rossokhata, and B. Bandyrskii, Functional-discrete method for an eigenvalue transmis- 
sion problem with periodic boundary conditions, Computational Methods in Applied Mathematics, 45 
(2005), no. 2, pp. 201-220. 

[16] M. Mikhailov and M. N. Ozisik, Unified analysis and solutions of heat and mass diffusion, Mir, John 
Wiley and Sons, 1984), 524R 

[17] M. Picone, Sur un probleme nouveau pour V equation lineire aux derivees partielles de la thebrie 
mathematique classique de I'elasticite, Colloque sur les equations aux derivees partielles, (Mai 1954). 

[18] N. Rossokhata, Analysis of an eigenvalue transmission problems with FD-method, Bulletin of the Uni- 
versity of Kiev, (2006), no. 1, pp. 194-203. 

[19] N. Rossokhata, FD-method for eigenvalue transmission problems with potential in space L 1; Bulletin of 
the University of Kiev, 3 (2007), pp. 161-168. 

[20] M. Schechter, A generalization of the problem of transmission, Annali della Scuola Normale Superiore 
di Pisa: Scienze Fisiche e matematiche, 14 (1960), no. 2, pp. 207-236. 

[21] V. Seng, K. Abbaoui, and Y. Cherruault, Adomian's polynomials for nonlinear operators, Math. Corn- 
put. Modelling, 24 (1996), no. 1, pp. 59-65. 

[22] G. Stampacchia, Su un problema relativoalle equazioni di tipo elliptico del secondo ordine, Richerche di 
Matematik, 5 (1956), pp. 3-24. 

[23] L. Vulkov, Well-posedness and monotone iterative method for a nonlinear interface problem on disjoint 
intervals, American Institute of Physics, Proceedings Series, 946 (2007). 

[24] P. Zhidkov, Basis properties of eigenf unctions of nonlinear Sturm- Liouville problems, Electronic Journal 
of Differential Equations, 2000 (2000), no. 28, pp. 1-13. 



